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Projected Entangled Pair States (PEPS) provide a framework for the construction of models where 
a single tensor gives rise to both Hamiltonian and ground state wavefunction on the same footing. A 
key problem is to characterize the behavior which emerges in the system in terms of the properties 
of the tensor, and thus of the Hamiltonian. In this paper, we consider PEPS models with Z 2 on¬ 
site symmetry and study the occurence of long-range order and spontaneous symmetry breaking. 
We show how long-range order is connected to a degeneracy in the spectrum of the PEPS transfer 
operator, and how the latter gives rise to spontaneous symmetry breaking under perturbations. We 
provide a succinct characterization of the symmetry broken states in terms of the PEPS tensor, 
and find that using the symmetry broken states we can derive a local entanglement Hamiltonian, 
thereby restoring locality of the entanglement Hamiltonian for all gapped phases. 


I. INTRODUCTION 

Correlated quantum many-body systems exhibit a 
wide range of unconventional phenomena. Their rich 
physics emerges from the intricate entanglement struc¬ 
ture of those systems, which however at the same time 
renders their theoretical study a challenging task. Some 
of the most important insights into the physics of those 
systems have thus been obtained through wavefunction 
ansatzes, such as the BCS or the Laughlin state. In re¬ 
cent years, ideas from quantum information have led to 
classes of ansatz wavefunctions constructed to capture 
the entanglement structure present in interacting quan¬ 
tum systems M- In particular, Projected Entangled 
Pair States (PEPS) [2j describe correlated many-body 
systems by associating a local tensor to each site which 
builds up the entanglement in the wavefunction through 
auxiliary indices. PEPS provide a faithful approxima¬ 
tion for low-energy states of systems with local interac¬ 
tions a a, making them the basis of powerful varia¬ 
tional algorithms @0. At the same time, they form a 
framework for the construction of models: From the local 
tensor, one can construct parent Hamiltonians which in¬ 
herit its symmetry structure and have a global wavefunc¬ 
tion built from the tensor as their ground state 00- 
Thus, they can be used to construct “PEPS models” 
(Fig. [lji), where the physics is encoded in a tensor A 
which gives rise to a local Hamiltonian H and a ground 
state wavefunction |?/>) on the same footing, forming a 
versatile framework for the study of correlated quantum 
systems 0-0] . Within this framework, the key problem 
is to understand how the local properties of the tensor 
(and thus the Hamiltonian) determine the global proper¬ 
ties, i.e., quantum order, of the wavefunction. 

In the last years, considerable progress has been made 
in the study of topological order 0] in PEPS models: 
It has been understood how topological order is related 
to the structure of the local tensor, and how this allows 
to construct all topological ground states from a single 
tensor 0 , 0 , 0 ] . More recently, the mechanism behind 
topological phase transitions within this framework has 
been clarified 0], and it has been found that it can be 


related to symmetry breaking in the so-called transfer op¬ 
erator 0], Finally, PEPS models have also been shown 
to provide a natural framework to study the entangle¬ 
ment properties of correlated quantum systems through 
entanglement Hamiltonians [271 ] associated to the bound¬ 
ary of the system 0 . 

The aim of this paper is to initiate the development of 
an analogous framework for studying conventional long- 
range order and spontaneous symmetry breaking within 
PEPS models, thereby complementing the understand¬ 
ing of topological phases. We consider PEPS models 
where the local tensor carries a Z 2 on-site symmetry, 
and investigate i) under which conditions long-range or¬ 
der emerges, and ii) how the different symmetry broken 
states can be described in terms of the single PEPS ten¬ 
sor. We first show that long-range order originates in 
a degeneracy of the transfer operator, and how the lat¬ 
ter in turn leads to spontaneous breaking of the symme¬ 
try under fluctuations. We subsequently characterize the 
symmetry broken states in terms of the PEPS tensor, 
and show that they correspond to the extreme points of 
the degenerate fixed point space of the transfer opera¬ 
tor. We demonstrate that the symmetry broken PEPS 
wavefunctions can be used to determine the spontaneous 
magnetization, and that they give rise to a local entangle¬ 
ment Hamiltonian, in contrast to previous findings [281 ] . 
which establishes that all gapped phases can be assigned 
a short-ranged entanglement Hamiltonian. 



FIG. 1. (a) PEPS models use a tensor A to provide a descrip¬ 

tion of a wavefunction \ip) and an associated parent Hamilto¬ 
nian H on the same footing, (b) Properties of the tensor such 
as on-site symmetries are inherited by H. (c) The transfer 
operator T encodes all correlations of the wavefunction. 
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II. PROJECTED ENTANGLED PAIR STATE 
MODELS 

Let us start by introducing PEPS. W.l.o.g., we restrict 
to a square lattice with Nh x N v sites. A PEPS model 
is described by a 5-index tensor A l a p i = 0,..., d — 1, 
a, /3, 7 , 5 = 0,..., D — 1, with d the physical dimension at 
each site and D the bond dimension. The PEPS wave- 
function I'f') = Y c ii...ijvl*i7 • ■ •) *iv} i s obtained by ar¬ 
ranging the tensors on the lattice and contracting the 
virtual indices a ,..., <5 of adjacent tensors as indicated 
by lines in Fig. [l}i, yielding We will consider 

systems on a long cylinder (or torus) of circumference 
N v . such that the dependence on the boundaries becomes 
negligible. 

To any PEPS, one can construct so-called parent 
Hamiltonians which have the corresponding PEPS as its 
ground state. Parent Hamiltonians ensure that the wave- 
function on a small (e.g., 2x2) patch looks correct, i.e., 
as being built from the tensor A, which is e.g. accom¬ 
plished by a projector onto the span of the tensors on 
that patch; under generic conditions, these Hamiltonians 
have a unique finite-volume ground state a By block¬ 
ing all N v tensors in a column, a PEPS can equally be 
considered as a quasi-ID tensor network (a Matrix Prod¬ 
uct State, MPS) with tensors B l and bond dimension 
D Nv . In MPS, a central role is played by the transfer op¬ 
erator T = Yi B 1 G B l \ in particular, appears in cor¬ 
relation functions at distance l, whose decay is therefore 
governed by the spectrum of T. For PEPS, T has itself a 
ID structure, cf. Fig. [lj:. This leads to a more complex 
behavior, since in the thermodynamic limit the system 
grows simultaneously in both directions, and thus, the 
dimension of the space on which T acts grows exponen¬ 
tially with the system size. 

Let us now turn towards systems with on-site sym¬ 
metries. In PEPS models, symmetries are encoded lo¬ 
cally in the tensor: A symmetry action u g on the physi¬ 
cal level translates to an action U g on the virtual sys¬ 
tem in a way where it cancels out when contracting 
tensors, see Fig. [lj:, givin g ri se to an invariant wave- 
function |T) = U® N |'L) [29(1 ■ This induces a sym¬ 
metry [T, Uf Nv ® Uf Nv ] = 0 of T, which is there¬ 
fore block-diagonal in a basis of irreducible represen¬ 
tations (irreps) of U g . At the same time, the parent 
Hamiltonian H enjoys by construction the same sym¬ 
metry, i.e., we obtain a PEPS model with symmetry 
[H 1 uf N ] = 0. An instructive example is the “Ising 
PEPS” [U with A = |O)(0,0,0,0| + |1)(0,0,0,0|, where 
the ket (bra) corresponds to the physical (virtual) indices, 
and |0) = cos0 |0) + sin0 |1), |0) = sin0 |0) + cos0 |1). 
This model has a Z2 symmetry with non-trivial action 
ui = U\ = X = (^ J). The wavefunction is of the form 

|*>= e _/3/2ff ° l(il ’'"’ ijv) |* 1 ,..., i N ) 

with H c 1 the classical 2D Ising model; it thus has the 
same er~ correlation functions and therefore undergoes a 


second-order phase transition at 


0 C = i arcsin 


1 

Vl+i/2- 


0.349596 . 


(1) 


III. LONG-RANGE ORDER AND 
SPONTANEOUS SYMMETRY BREAKING 

As we have seen, PEPS form a natural framework to 
model systems with on-site symmetries: By encoding 
the symmetry locally into the tensor, we obtain a lo¬ 
cal Hamiltonian with the same symmetry which depends 
smoothly on the parameters of the tensor, and whose 
ground state wavefunction can be constructed from the 
very same tensor. However, we have also found that by 
changing parameters, such a model can undergo a phase 
transition to a symmetry broken phase. Nevertheless, 
the PEPS wavefunction |'L), which is the unique ground 
state of the system, remains invariant under the symme¬ 
try throughout the phase diagram, |T) = it® JV |4'); e.g., 
for the Ising PEPS at zero temperature (0 = 0), we have 
I’L) = -^-(|0• • • 0) + |1 ■ • • 1)). This leads to the central 
question of this work: How can we characterize symmetry 
breaking in PEPS models, and how can we construct the 
distinct symmetry broken states starting from a single 
symmetric tensor A? 

Let us for a moment leave aside PEPS models and 
consider the analogous question for general Hamiltoni¬ 
ans with a symmetry: Generically, these systems have a 
unique finite volume ground state which therefore carries 
the full symmetry of the Hamiltonian. Since the average 
magnetization (jj Yi %i) i n those states is zero, the or¬ 
dered phase is characterized through non-vanishing long- 
range order (-Tj Yij Zj ) ~ const. In order to obtain 
the symmetry broken states, one couples the system to 
a small external field H' = H + h Y-, and considers 
the ground state , taking first the limit N —> 00 

and subsequently h —> 0; if in the limit ; ^ Yi %i) ^ 0, 
one says that the system exhibits spontaneous symmetry 
breaking. 

How are these two notions of long-range order and 
spontaneous symmetry breaking related, and can one 
construct an approximate symmetry broken state for a 
finite system? Given a unique finite-volume ground state 
with long-range order such as |’L) « |0 • • • 0) + |1 • • • 1) 
and order parameter Z, one can construct an orthogo¬ 
nal state |4>) oc Yi^iY^) whose energy approaches the 
ground state as IV —> 00 [3j]]. A generic local pertur¬ 
bation will select a unique ground state from the space 
spanned by |\I/} and |<f>), assuming a gap above. In fact, 
it can be proven that the only states which do not exhibit 
long-range order for any observable O with (O) = 0, i.e., 
which are stable under arbitrary perturbations and thus 
form the symmetry broken states, are |\&) ± |$) [32j]; this 
is, the symmetry broken states can be constucted from 
the symmetric ground state |'L) alone. It is straightfor¬ 
ward to see from this argument that the value of the 
long-range order is just the spontaneous magnetization 
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squared. It should be pointed out, however, that the 
above argument relies on the validity of the perturba¬ 
tive treatment which cannot be justified rigorously in the 
limit N —> oo, as the total perturbation diverges, and in 
fact, the equality between the two notions of symmetry 
breaking and the values of the corresponding order pa¬ 
rameters has not been proven rigorously except for very 
few cases such as the Ising model |33|. 

In the following, we will follow a very similar reason¬ 
ing in order to relate long-range order and spontaneous 
symmetry breaking in PEPS, and in particular to show 
how to construct the set of symmetry-broken states (i.e., 
those which are stable under arbitrary local perturba¬ 
tions) starting from a single object, namely the unique 
symmetric PEPS wavefunction |\H) and its associated 
symmetric tensor A. Here, a central role will be played 
by the one-dimensional transfer operator and in particu¬ 
lar its spectral properties and eigenstates, partly taking 
the role of the Hamiltonian in the preceding discussion. 
We will start by showing that long-range order in the 
wavefunction is closely related to an approximate degen¬ 
eracy in the spectrum of the transfer operator. Subse¬ 
quently, we study the behavior of the transfer operator 
under perturbations, and we show that arbitrary pertur¬ 
bations induce a splitting of the degenerate subspace in 
a fixed, basis- independent of the perturbation —which 
therefore describes the symmetry broken states of the 
model. 

In the following discussion, we will restrict to the sym¬ 
metry group Z 2 , and denote the physical (virtual) sym¬ 
metry action by ir (X). 


assume this in the following. (Corresponding restrictions 
on Nh will not be necessary in the limit Nh —> 00 away 
from fixed point wavefunctions.) 

In ID, i.e., MPS, long-range order is in one-to-one cor¬ 
respondence with a degeneracy of the transfer operator 
T—roughly speaking, the only way to build long-range 
correlations over arbitrary distances using the constant 
virtual dimension is an exact degeneracy of the lead¬ 
ing eigenvalue of T. In 2D, however, the situation is 
much less clear due to the exponentially growing dimen¬ 
sion of the virtual space on which T acts. However, as 
we will show, under certain conditions Eq. fl5j) still im¬ 
plies the existence of an almost degenerate second eigen¬ 
value. To this end, define a “dressed” transfer operator 
T z„, where we insert Z on the physical index at position 
s = 1 ,...,N V , and correspondingly ^ X)«=i 5 

such operators arise in expectation values such as Eq. ©• 
Assume that T is diagonalizable, T = ^)Ai| r i){h\ with 
(k\rj) = Sij. Here, we use |-) to denote (eigen-)vectors 
on the level of the transfer operator, i.e., on the virtual 
indices; due to the ket-bra structure of the transfer op¬ 
erator, these vectors can themselves be regarded as op¬ 
erators, in which case we will write them without brack¬ 
ets. Since [T, <g) X® Nv ] = 0, the eigenvectors of 

T transform even or odd under the symmetry (i.e., as 
different irreps), and we will denote the largest eigen¬ 
value in the even (odd) symmetry sector by A+ (A_), 
with corresponding eigenvectors |?"±)(Z±|. We thus have 
that (X®^” ®X'® JV ”)|r±) = ±|7-±), or (if interpreting the 
eigenvectors as operators) 


A. Long-range order and the spectrum of the 
transfer operator 

Let us first see how long-range order in a PEPS model 
with Z 2 on-site symmetry is related to the properties of 
the underlying tensor and in particular the transfer oper¬ 
ator. We start by defining long-range order (in analogy 
to Ref. [3^|i : We say that |^) has long-range order if 
there exists a local operator Z = with Zx = —xZ 
and ||Z||op < 1 (the order parameter ), and c > 0 s.th. for 
sufficiently large N v , 

cK (2) 

u m,n 

where the sum runs over all sites—i.e., on an infinite 
cylinder, the spins are correlated at least over a dis¬ 
tance proportional to the circumference. [Away from 
fixed point wavefunctions, we indeed cannot expect cor¬ 
relations along the cylinder over arbitrary distances; on 
the other hand, e.g. in the Ising PEPS correlations only 
break down after a distance exp(ciV„).] Note that Z can 
depend on the lattice site, such as in the case of an anti- 
ferromagnet. In such a case, we will also have to restrict 
the possible values of N v accordingly, e.g. to even num¬ 
bers (or we have to block tensors), and we will tacitly 


X® jV “r±X® JV ” = ±r± . (3) 

Note that this implies that r_ cannot be positive semidef- 
inite, and since T must have a positive semidefinite fixed 
point (as it is a completely positive map), it follows that 
|A+|>|A_|. 

Our goal is to show that A+ and A_, as a function 
of N v , become degenerate as N v —y 00 . We can thus 
assume that |A_| < |A+| (otherwise there is nothing to 
show). We will additionally restrict to the case where the 
largest eigenvalue A + is non-degenerate also within the 
even parity sector, since an exact degeneracy at finite N v 
hints an additional symmetry in the wavefunction which 
we would have to incorporate for a full description. This 
implies A+ > 0 (HI, [il5j and we can w.l.o.g. normalize our 
tensors such that A+ = 1. We can now re-express Eq. © 
as 


tr\T 7 TPTrT Nh -P- 2 } 

cN v < lim Y, 1 + 1 > 


Nh—>oo 


p —0 


tr[T JVfc ] 


(4) 


where the additive term +1 stems from the case where 
both Z 's are in the same column, and using ||Z|| op < 1. 
In order to simplify the r.h.s., we split the sum at p = 
Nh /2 — 1 and use cyclicity of the trace, together with 
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p £->■ Nh — p — 2 , to obtain two identical sums, 


B. Spontaneous symmetry breaking 


tr T yT p TyT Nh ~ p ~ 2 \ 

cN v < lim 2 > —^-- + 1 . 

- jv„->oc ^ tr 

p— 0 L J 


(If Nh is even, the term p = Nh/ 2—1 appears only once, 
but this bears no relevance.) Since the largest eigenvalue 
of T is non-degenerate, we have that 

||T M -|r + )(4||| tr <cr M , (5) 


where F < 1 upper bounds the second largest eigenvalue 
of T [Sfij]. (Note though that c can heavily depend on 
N v and properties of T.) Using a sequence of triangle 
inequalities and taking the limit Nh —> oo (see Appendix 
for details), we can then replace T Nh —> |r+)(4l and 
obtain 


cN v < 2 £tr[T^TH4|r + )(4l]+l ( 6 ) 

p—0 

oo 

= 2j2J2 X ^ l +\ T z\rmTz\r + )+l , 

- - v „ ✓ 

i p —o 

=:rrii 

where in the last step, we have expanded T in its eigen- 
basis. 

As anticommutes with the symmetry, (Z + |T^|rj) = 
0 for eigenvectors | rt) from the even sector. Thus, only 
eigenvectors from the odd sector contribute to the sum. 
For those eigenvectors, |Aj| < |A_|, and therefore 


p 

If now nrii > 0, we have = z\ r +) — 1 

(since ||Z||o P < 1 ) and thus 


< 


l-|Ai 


rad < 


1 — I A—| 


CNy < 


2 

iHAii 


+ 1 ) 


i.e., |A_ | > 1 — 0(1/N V ): We find that long-range order 
implies that the gap of T closes in the thermodynamic 
limit at least as 1 /N v . The required condition to 7 > 
0 is automatically satisfied if T and T g are hermitian, 
which in particular holds if the tensor and thus the PEPS 
model is invariant under combined reflection and time 
reversal, but turns out to be true also for a range of 
other models [37| . 

We thus find that under certain conditions, long-range 
order implies that the transfer operator has approxi¬ 
mately degenerate eigenvalues A+ « A_ in the even and 
odd sector, together with non-vanishing matrix elements 
(4|T^|rrp). Conversely, an approximately degenerate 
eigenvalue in the odd sector, together with non-vanishing 
matrix elements, will clearly give rise to long-range order. 


Let us now investigate spontaneous symmetry breaking 
in PEPS models with long-range order, and in particular 
how to construct the symmetry broken states. The sym¬ 
metry broken states are those states which remain ground 
states under generic perturbations, where for a given per¬ 
turbation, we first need to take the thermodynamic limit 
and subsequently take the strength of the perturbation to 
zero. In the context of PEPS models, the natural pertur¬ 
bations to consider are perturbations of the tensor which 
can be realized by acting solely on the physical index, 
AapryS > XX ^ + ^)ijA j a p lS with || A|| < 1 , as these both 
correspond to small perturbations of the parent Hamil¬ 
tonian and leave us in the PEPS manifold [38j. Such a 
perturbation gives in turn rise to a perturbation of the 
transfer operator, T —> and thus the symmetry bro¬ 

ken states will exactly correspond to boundary conditions 
given by the stable fixed points of T[aj in the appropriate 
limit. 


1. Hermitian transfer operator 

How will the transfer operator of a system with long- 
range order respond to such perturbations? To this end, 
consider the limit N v —> oo where A_/A+ —> 1, and let 
us assume that T has a gap below the two degenerate 

eigenvalues, T = |r+)(4| + l r -)(M +_A sufficiently 

small perturbation will then induce a splitting within this 
subspace without mixing it with lower-lying eigenvectors; 
the basis in which this splitting occurs is independent of 
the strength of A and thus leaves us with a single fixed 
point T^j- ~ li m it Nh —> oo, even if later 

|| A|| —> 0. This fixed point is therefore the boundary con¬ 
dition for the symmetry broken state under perturbation 
A. 

Naively, one might expect that the symmetry broken 
fixed point sensitively depends on the exact form of the 
perturbation chosen. However, as we will show in the 
following, the symmetry broken fixed point 4'!')(41 ^ n_ 

dependent of A. We will first consider the case of a her¬ 
mitian transfer operator, i.e., (r 7 | = (41 and ( 4 r j) = ^iji 
and a hermiticity-preserving perturbation. A generic 
such perturbation to T will induce a splitting of the two 
degenerate eigenvalues in an orthogonal basis 

kt) « l r +) +7k-) 

K) oc |r+) - —|r_) . 

Since are non-degenerate eigenvectors of T[a], which 
is a completely positive (CP) map, they must be hermi¬ 
tian, and thus 7 £ R. Since T[aj is a completely positive 
map, its leading (non-degenerate) eigenvector must be a 
positive operator, 


r-f. oc r + + qr_ > 0 . 


(7) 
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Using Eq. © and the fact that conjugation preserves 
positivity, we find that also 

r+ - qr_ = X 0JV ” (r + + > 0 . (8) 

Since for a pair of positive operators P, Q > 0 , it holds 
that tr[PQ] > 0, it follows that 

0 < tr[(r+ + yr_)(r+ — 7r_)] 

= ir+V+) - l 2 {r-\ r ~) , 

and thus |7| < 1 . By changing the sign of A, we can 
exchange the two leading eigenvectors, i.e., I77) becomes 
the leading eigenvector, and thus 77 > 0 . Following the 
same line of reasoning as before, this yields |y| > 1, and 
thus, 7 = ± 1 . We thus find that regardless of the pertur¬ 
bation, the transfer operator always acquires the same 
pair of fixed points 

kn) c* k+) ± k-), 

solely as a consequence of its complete positivity. Note 
that 7"i7 are the extremal positive states r + + yr_ > 0, 
i.e., the stable symmetry broken states are those where 
the symmetry is maximally broken, as intuitively ex¬ 
pected. 

2. Non-hermitian transfer operator 

Let us now consider the non-hermitian case. Since 
k+)(4l + k-)(M = ktXkl + ktXkl, it; follows that 

ku) oc k+) ± A k-) and (hi I oc (41 ± ik-l 

where hermiticity of the eigenvectors implies Asl, and 
we w.l.o.g. choose A > 0 . Again, complete positivity 
of T[aj and T '4 , together with the possibility to change 
the ordering of eigenvectors by changing the sign of A, 
implies that 774., Z47 > 0 . To determine A, consider 

S := suppr + = (kerr + ) J - , ( 9 ) 

and let II5 be the orthogonal projector onto S. Com¬ 
plete positivity of T implies that supp r_ C S f39j , and 
therefore supp 774. C S , i.e., ILsr-iqlLs = 774,. Define 

hi = ns4i4s > 0 ; 

since (ktktt) = = tr[ n 5 Z^ n 5r t 4.] = 

tr^rpiJ = (l n \r n ) = 1, we have that 

hi ± 0 • (! 0 ) 

As 774. = r+ ± Ar_, and r_j_ > 0 , A > 0 , it follows that 
ker 77 fi S (ker 77 (~l S') is contained in the negative (posi¬ 
tive) eigenspace of r_, which implies that 

(ID 


Moreover, since tr^r^] = tr^j^t] = (ktktt) = we 
have suppZ^ C kei' r^-. and thus 

supp hi C ker rjq n S (12) 

(as suppZ -17 C S). Eq. has two implications: First, 
together with Eq. m,hi ^ 0 , it shows that ker ^05 ^ 
0. Since S = supprp [Eq. ©], this implies that 777 are 
again the extremal points of the positive cone rp + Ar_ > 
0. Second, by combining Eq. m with Eq. (ED, we find 
that supp/-)- _L supp/ 7 , and thus 

0 = tr[Z7Z7] = tr[II 5 (4 + AL)n s (4 - XI-)] , 

which allows us to determine the corresponding value 

of A as A = [tr[(ILs4) 2 ]/tr[(II < 5 Z_) 2 ])] Together, 
this proves that also in the non-hermitian case, the 
fixed points of the transfer operator under perturbations, 
i.e., the symmetry broken states, are uniquely deter¬ 
mined independent of A, and correspond to the extremal 
symmetry-broken states. 

IV. NUMERICAL STUDY AND 
ENTANGLEMENT HAMILTONIANS 

Having understood the structure of the symmetry bro¬ 
ken states in the thermodynamic limit, let us now inves¬ 
tigate the accuracy of our result for finite systems, and 
use our findings to re-examine the entanglement spectra 
and Hamiltonians of models with symmetry breaking. 

To this end, we have performed numerical simulations 
for two models. First, we have considered the Ising PEPS 
introduced in Sec. El As its transfer operator is isometric 
to the transfer operator of the 2D classical Ising model 
and it exhibits the same a z correlation functions, it allows 
us to benchmark our numerical findings against exact 
results. 

As a second model, we have studied the square lat¬ 
tice AKLT model with nematic field |28[- The AKLT 
model HI is constructed by placing spin- d singlets on 
the links and subsequently projecting them onto the spin 
2 subspace, i.e., A is of the form ns =2 (l Cg> 1 ® a y ® a y ), 
where ns =2 projects onto the spin -2 space; the model 
can be mapped to a ferromagnetic model through a sub¬ 
lattice rotation. Subsequently, a “nematic field” A —> 
X) NijA^p 7<5 , N = exp(aS'4 favoring large values of 
is applied to the AKLT tensor which yields a model with 
1 /( 1 ) xi Z 2 symmetry and eventually leads to breaking of 
the Z 2 symmetry [28| . 

In all cases, the simulations have been carried out 
through an exact diagonalization of the transfer oper¬ 
ator (cf., e.g., Ref. [40|), where the transfer operator T 
is applied to a vector by sequential contraction of indices 
rather than by building its matrix representation. Thus, 
the largest objects which need to be stored are on the 
order of the size of an eigenvector, D 2Nv , rather than 
that of the transfer operator, £) 4Ar "; at the same time, 


(ker qn 5 ) T (ker 77 fl S) . 
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FIG. 2. Numerical results for the Ising PEPS. (a,b) Er¬ 
ror analysis for A = tZ, 9 = n/16: a) Scaling of the error 
(5 = 1 — \(f\g)\ 2 between the fixed point of T[a], |/), and 
of its projection onto |r±), | g)\ we find 5/N v ~ e 2 . (The 
data for different N v almost exactly coincide.) b) The devi¬ 
ation of u) = |(r_ \g)/(r+\g)\ 2 from 1 vanishes exponentially 
in N v . (c) Magnetization, computed from the symmetry bro¬ 
ken states |r+) + |r_), for N v = 4,. .., 12 (blue), extrapolated 
curve (red), and exact solution from the correspondence to 
the classical 2D Ising model (dashed green). 

one obtains results with machine precision, allowing for 
an accurate scaling. 

A. Error scaling 

Our derivation involved a number of assumptions, in 
particular neglecting the coupling to eigenstates other 
than |r±) in second and higher order perturbation theory 
even in the limit N v —> oo, and neglecting the splitting 
between A± [4lt | . In order to assess the validity of neglect¬ 
ing the coupling to the other levels, we have compared the 
fixed point |/) of the perturbed transfer operator IW, 
with A a general perturbation, to the fixed point | g) ob¬ 
tained after projecting T[aj onto |r±) and (l±\. Fig. [2b 
shows the result for the Ising PEPS, where we find that 
the error S = 1 — |(/|ff)| 2 scales as <5 ~ -/V„||A|| 2 . This scal¬ 
ing is consistent with a second order perturbation treat¬ 
ment, where one assumes that the local perturbation A 
has only a local effect on |r±), and thus can only be un¬ 
done by a term in its vicinity, leading to 0(N V ) contribu¬ 
tions, rather than 0(7V 2 ), of amplitude ||A|| 2 each. The 
same behavior is found for the nematical AKLT model, 
cf. Fig. [3b. We have subsequently also verified the valid¬ 
ity of neglecting the splitting between A±: Fig. [5)3 shows 
that the exact solution for the two-dimensional degener¬ 
ate subspace, | g), in turn converges exponentially to |rqq) 
as N v —> oo, consistent with an exponentially vanishing 
splitting 1 — A_/A+. 

B. Magnetization 

The knowledge of the boundary conditions |r-|q) and 
(l- tq| of the symmetry broken states allows us to com¬ 



FIG. 3. Nematic AKLT model. (a) Error analysis, 
cf. Fig. [2jr. (b) Critical exponents for the correlation length, 
determined as £ gap = —1/ln(A_/A+) in the disordered phase 
and £mag = £(0) from the magnetization fit (cf. Fig. [2b) 
in the symmetry-broken phase, (c) Magnetization m for 
N v = 4,..., 12 and extrapolated value (cf. Fig. [2b); the inset 
shows the critical exponent of m. 


pute the approximate spontaneous magnetization for fi¬ 
nite N v . For the Ising PEPS, we find that the finite N v 
data reproduces the analytic expression for the infinite 
system very well away from the critical point, with the 
error vanishing exponentially in both N v and the distance 
from the critical point; this is consistent with the expec¬ 
tation that the error should vanish as e~ Nv ^ with £ the 
correlation length. On the other hand, there are strong fi¬ 
nite size effects as we approach the critial point, as shown 
in Fig. [2b; indeed, it is natural to expect that finite size 
effects will dominate as soon as £ ~ N v . Nevertheless, 
it is still possible to accurately extrapolate the magneti¬ 
zation curve to the thermodynamic limit even very close 
to the critical point. To this end, we first determine the 
behavior of the correlation length in the proximity of the 
critical point by fitting mg(N v ) = a + bexp(—N v /£g) as 
a function of N v , and subsequently fitting l/£e quadrat- 
ically (shown in the inset; note that this makes an im¬ 
plicit assumption about the critical exponent of £); the 
critical point from the fit is 0 C ^ ps 0.349666, as com¬ 
pared to the exact value 0 Ci6X ps 0.349596. The resulting 
fitting function £'(0) is then used in the scaling ansatz 
mg(N v ) = me{oo) + fN~ 9 exp[— N v /£'(0)], where the al¬ 
gebraic factor N~ 9 accounts for short-range effects [42j |. 
Altogether, this extrapolation yields a very accurate ap¬ 
proximation to the analytical curve even in the regime 
very close to the critical point where finite size effects 
seemingly dominate, as shown in Fig. [2b. 

We have subsequently applied the same analysis to de¬ 
termine the magnetization curve of the nematic AKLT 
model. The data for N v = 4,..., 12 and the extrapo¬ 
lated magnetization curve are shown in Fig. [3bj and we 
find a value of a c ps 0.0447 for the critical point. Us¬ 
ing this value of a c , we have subsequently studied the 
critical scaling of the model. Fig. [3b» shows the criti¬ 
cal scaling of two different correlation lengths—£ m ag, ob- 
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FIG. 4. Interaction strength vs. interaction range of 
H ent [28} for (a) the Ising PEPS (9 = 0.19, 0.23,..., 0.51) 
and (b) the nematic AKLT model (a = 0,0.01,..., 0.1), 
where blue=trivial phase, green=symmetry broken phase, 
and red=critical point. The insets show H en t computed from 
<7+ , where the interaction becomes long-ranged after the phase 
transition. The main panels show H en t derived from the sym¬ 
metry broken states (gray lines indicate the data from <r+), 
and we find that in both cases, H ent becomes more local again 
after the phase transition. 


tained by fitting the magnetization as a function of N v 
as m(N v ) = a + 6exp(—-/V„/£ mag ), and £ gap , obtained 
by an exponential fit £{N V ) = £ gap + a exp(—bN v ) of the 
correlation length £(N V ) = — ln(A_/A+) extracted from 
the gap of the transfer operator—both of which yield a 
critical exponent of v ~ 1 for the correlation length. The 
inset of inset of Fig. [3]c shows the corresponding analysis 
for the magnetization, which is compatible with a critical 
exponent /3 = 1/8, therefore suggesting that the nematic 
AKLT model is in the 2D Ising universality class. 


C. Entanglement Hamiltonian 


Let us now turn towards entanglement spectra and 
the entanglement Hamiltonian. The entanglement spec¬ 
trum of a state on e.g. an infinite cylinder is given by 
the spectrum of the reduced density operator pl of the 
half-infinite cylinder, and its low-energy part can be in¬ 
terpreted as the spectrum of a Gibbs state e~ H<,nt of an 
“entanglement Hamiltonian” H ent [271] . In the context of 
PEPS models, it has been shown (28| that the spectrum 
of pl is equal to the spectrum of the symmetrized fixed 


point of T, 



where • denotes a possible label of the fixed point, such as 
in the case of multiple sectors. In turn, er, allows to define 
the entanglement Hamiltonian through H ent := — In er,. 

In Ref. 0, H ent has been studied for both the Ising 
PEPS and the nematic AKLT model based on the unique 
fixed point of T for finite N v , i.e., • = +, and it has been 
found that it is (quasi-)local (i.e., the interactions de¬ 
cay exponentially with distance) in the trivial phase, the 
interaction length diverges at the phase transition, and 
H ent remains long-ranged in the symmetry-broken phase; 
the corresponding results for both models are shown in 
the insets of Figs. QJi,b. 

However, as we have argued, in the symmetry broken 
phase the physically meaningful states, i.e., those which 
are stable under arbitrary perturbations, are the sym¬ 
metry broken ones, and we should therefore rather use 
the symmetry broken fixed points • =t,4- when deter¬ 
mining iTent. In Fig. [H we compare the interaction range 
of Hent derived from • =t,4- [43 with the one obtained 
from • = + (insets) for the Ising PEPS (Fig. |4ji) (44| 
and the nematic AKLT model (Fig. 0]:>): We find that 
by considering the latter, the locality of H ent in the sym¬ 
metry broken phase is restored, and H ent diverges only 
at the phase transition, in accordance with the intuition 
that the interaction length of the entanglement Hamilto¬ 
nian should reflect the characteristic length scale of the 
system, and thus should be finite away from the critical 
point. 

Together with previous findings on how to recover lo¬ 
cality of H e nt in topological phases [25[, this shows that 
by considering the physically meaningful (this is, stable) 
fixed points, a local entanglement Hamiltonian can be 
obtained for all gapped phases, as intuitively expected. 


V. CONCLUSIONS 


We have studied the occurence of long-range order and 
spontaneous symmetry breaking in PEPS models. We 
have shown that long-range order is closely related to 
a degeneracy in the transfer operator, and have char¬ 
acterized the symmetry broken (i.e., stable) states in 
terms of the fixed points r± of the transfer operator, 
which we found to be the extremal positive semi-definite 
states r-jq = r + ± Ar_ > 0 irrespective of the model 
and the perturbation. These fixed points do not only 
yield the physically relevant boundary conditions in the 
symmetry broken phase, but each one by itself already 
carries the full information about both fixed points, as 
r± = r-fzLX® N v r -{X® Nv . Moreover, they give rise to local 
entanglement Hamiltonians, thereby establishing locality 
of the entanglement Hamiltonian for all gapped phases. 
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vector of r+ + r_, but we have used |r) = |1... 1) which 
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Appendix: Relation between long-range order and gap of the transfer operator 


In this appendix, we give a detailed derivation of the equality of the r.h.s. of Eqs. and ©. We start by splitting 


N h -2 

E 


p =o 


tv[T zTPJ zY Nh ~P- 2 ] 
tr [T Nh ] 


5(0, N cut ) + S(N cut + 1, N h — 2) , 


where N cut = L^rJ ~ 1> and 


S(a, b) 


E 

p=a 


tr[T^T p T^T Ar '* _p_2 ] 

Hffwj 


Due to cyclicity of the trace, S(N cut + 1, Nh — 2) = 5(0, N cut + n), where N cut + k = Nh — 2 — (7V cut + 1) and thus 
k = —1, 0, depending whether Nh is even or odd. 

We will now show that 


v lim^5(0,7V C ut + «) = E tr [T/PT*|r + )(4|] ■ 

p—0 


To this end, we will make use of Eq. ©, 

l|T M — |r_|_)(Z + | || tr < cT m © 

where T < 1, as well as 


| tr[T^T^]| < ||T^|| op ||T p || op ||T^|| op ||X|| tr < C||X|| tr (A.l) 

with £ := (c+ l)||T^|| 2 p , which can be shown using Holder’s inequality, the submultiplicativity of the operator norm, 
and 

© 

IITlop < ||T P - |r + )(4l|| op + |||r + )(4l|| op < ||T p - |r+)(4l|| tr + 1 < + 1 < c + 1 , 

and finally 

| tr[T JV '*]| > |tr[|r+)(4|]| - | tr[T^ - |r+)(4|]| > 1 - \\T n * - |r+)(4| || tr > 1 - cT^ . (A.2) 

We now have 


tr[T^TPT^T Nh ~P~ 2 ] 

-tr[T^'PT^|r + )(4|] 




tr[T^] 




tr[T ^T P T ^f Nh ~ p ~ 2 ] 

tr[T^TPT^|r+)(4|] 


tr pI4T p H41 f +) (41 ] 

tr[T^TPT^|r + )(4|] 

trJT^] 

tr[T^] 


tr[T N »] 

tr [| r +)(4|] 


[AIMM ( \\T Nh ~ p ~ 2 - |r + )(4||[ tr tr[|r + )(4|]-tr[T^] 

1-cT** +S||I +A t +l||tr tr[TW] t r[|r+)(4|] 


151X21 

< 


Ccrw-p-2 C T Nh 

l_ cr w fc + C||k+)(4l|| t r l _ cT N h 


< 2 (cvT Nh - p - 2 
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where in the last step we have assumed that Nh is sufficiently large such that 1 — cT Nh > and have introduced 
v := 1 + |||r+)(Z+||| tr . It follows that 


N c ut+zt 


5(0,JV cut + K)- £ tr [T^T^|r + )(Z + |] 


p—0 


< 


N c ut+zt 

E 

p—0 


< ^ x 2 CcvT Nh/2 


-1 


where we have used N cut + k < — 1 and Nh — p — 2 > Nh/2 — 1. Clearly, the rdr.s. goes to zero as Nh —> oo, and 

thus, 


lim 5(0, N cut + «) = [T^T^|r+)(Z + |] 

JV h ->oo —' 

p=0 


as claimed. 




